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Abstract. Although implicit-expHcit (IMEX) methods for approximating so- 
lutions to semilinear parabolic equations are relatively standard, most recent 
works examine the case of a fully discretized model. We show that by dis- 
cretizing time only, one can obtain an elementary convergence result for an 
implicit-explicit method. This convergence result is strong enough to imply 
existence and uniqueness of solutions to a class of semilinear parabolic equa- 
tions. 



1. Introduction 

The use of implicit-explicit (IMEX) methods for approximating semilinear par- 
abolic equations is well-established [1] . Many of the recent works on these methods 
employ discretizations in both space and time. These fully discrete approximations 
can be computed directly by a computer. However, one can obtain a stronger con- 
dition for convergence of the approximation if only the time dimension is discretized 
[2] . We show how an even stronger condition for convergence is met by the Cauchy 
problem for 



oo 



(1) ^^^^^^Au{x,t) + )^a,{x)u\x,t), 

where £ L^(R") n L°°(R"'), and how convergence of this method provides an 
elementary proof of existence and uniqueness of solutions. Existence and uniqueness 
of solutions for ([1]) under reasonable initial conditions have been known for some 
time. For instance, [3] and fB] contain straightforward proofs using semigroup 
methods. The purpose of this paper is to show how a more elementary proof can 
be obtained from a sequence of explicitly computed discrete-time approximations. 

The Cauchy problem for ^ arises in a variety of settings. Notably, some 
reaction-diffusion equations are of this form [3]. Another application is the spe- 
cial case 

fill ( If f\ 

— Au{x, t) — u {x, t) -\- ao(x). 



dt 

where ao is a nonzero ftmction of x. This situation corresponds to a spatially- 
dependent logistic equation with a diffusion term, which can be thought of as a toy 
model of population growth with migration. 
Following ^J, the approximation to be used is 



oo 



(2) M„+i = (/-/iA) ^(u„ -I- /i^a.O, 

which is obtained by inverting the linear portion of a discrete version of ([T]). For 
brevity, we shall call ^ the implicit-explicit method. (In the summary paper 
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[T], this is called an SBDF method, to distinguish it from other implicit-explicit 
methods.) One can compute the operator (/ — hA)~^ explicitly using Fourier 
transform methods, and obtain a proof of the numerical stability of the iteration 
as a whole. 



2. A VERSION OF THE FUNDAMENTAL INEQUALITY 

In order to simplify the algebraic expressions, we make the following definitions. 
Definition 1. Let 

oo 

(3) Fiuix,t)) = Auix,t)+Y, 

ai{x)u\x, t), 

i=0 

and 

oo 

(4) G(u(x,t)) = ^a,(x)u'(a;,<). 

4=0 

Definition 2. Define the analytic functions 

oo 

(5) gi{z) =^\\a^\\lz\ 
and 



(6) 500(2:) = XI 11°' I 



i=0 

Since we do not discretize the spatial dimension, we can employ some of the 
theory of ordinary differential equations. We therefore first prove a variant of the 
fundamental inequality for ([1]) as is done in [5j . The fundamental inequality gives a 
sufficient condition for approximate solutions to converge. A slightly weaker version 
of Lemma [3] was obtained in Theorem 3.1 of [2], where the existence of solutions 
was required. 

Lemma 3. Suppose {ui}°li is a sequence of piecewise functions Ui : [0,T] 
C2(M") n ii(R") n i°°(R"), such that 

(1) there exist A, B > so that for each i and t £ [0, T], |lui(t)||i < A and 

\\u^{t)\\oo < B, 

(2) for each i andtCz [^^T], the series gi{\\ui(t)\\i) and gac{\\ui{t)\\oa) converge, 

(3) for each t £ [0,T], ||^Mi(i) — F(ui{t))\\oo < o,nd limi^oo e-i = 0, and 

(4) ui{0) = UiiO) for alli>0 

Then for each t G [0,T], {ui{t)}°li is a Cauchy sequence in i^(R"). 

Proof. Let i,j > be given. Let rj{t) = \\ui{t) — Uj{t)\\2 — J {ui{t) ~ Uj{t))'^dx. 
Notice that the fourth condition in the hypothesis gives 77(0) = 0. 

rj'it) = 2 J (u'S) - u,{t)) Mt) - u,{t)) dx. 

But, ll^Ui(t)— F(iii(t))||oo < Ei is equivalent to the statement that for each i S [0,T] 
and X e M", 

F{ui{x, t)) - ei < u-(a;, t) < F{ui{x, t)) + e^. 
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giving 



V'{t) < j {F{ui{t))- F{uj{t))){ui{t)-Uj{t))dx + 2{€i + ej) j \ui{t) - Uj{t)\dx 

< (Aui(i) + G{ui{t)) - Auj{t) - G{uj{t))) {ui{t) - Uj{t))dx 
+2{ei + ej)\\ui{t)-Uj{t)\\i 

< 2 y" {/\{u^{t) - Uj{t))) {u^{t) - Uj{t))dx 

+2 j {Giuiit)) - G{uj{t))) {Ui{t) - Uj{t))dx + 2{ei + ej)\\ui{t) - Uj{t)\\i 

< -2 y" II V(wi(t) - Uj{t))fdx + 2\\G{uiit)) - G{uj{t))\\2\\ui{t) - Uj{t)\\2 

+2{ei + ej)\\u,{t) ~ Uj{t)\\i 

< 2\\G{ui{t)) - G{ujm\2\\ui{t) - n,(i)||2 + 2(e, + e,)||n,(t) - uj{t)\U. 

Now also 

\\G{ui{t)) - G{ujm\2 = 



fe=0 



fe=0 



< Ell«fc||ooW / {u't{x,t)-U^^{x,t)fdx 
fe=0 V 



< 



< 



fe=0 



'fc-1 



{u^ix,t)-ujix,t)f J2 uT{x,t)uY'''-\x,t) dx 



El 

fe=0 



fc-l 



m=0 



\\Ui{t) - u,{t)h 



< \Y,\\ak\\^kB^-^\\\ui{t) - Uj{t)\\2 



\k=0 



< g'M\\ui{t)-uM\2, 



which allows 



ri'it) < 2g'^{B)\\u,{t)-u,it)\\l + 2{ei + ej)\\u,{t)~u,it)\\,. 
< '^gUBHt) + 2(ei + ej)\\ui{t) - Uj{t)\U. 



r,'{t)-2g'^{B)r,{t) < 2{ei + ej)\\ui{t) - Ujit)\h 

I (r;(t)e-2«^(^)*) < 2(e, + e,)e-2^~(^)*||u.(t) - (t)|| i, 
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SO (recall 77(0) = 0) 

v{t) < 



2(e. + e,) / e-2ff:.(s)-||u,(s)-u^.(s)|jids 



< 



2{ei + ej) \\ut{s) ~ Uj{s)\\ids 



Hence as i, j — 
sequence in 



> 00, r]{t) for each t. Thus for each t, {ui{t)}°^^ is a Cauchy 
R"). □ 

Remark 4. Since C2(M")nLHR")ni°°(K") C L'^{W) and is complete, Lemma 
[3] gives conditions for existence and uniqueness of a short-time solution to ([1]). 

Lemma 5. Suppose {ui{t)}°^-^ is the sequence of functions defined in Lemma\^ 
and that u{t) — lim^^oo Ui{t) in L^(IR"). Then 

(7) u'{t,x) — lim u[{t,x) for almost every x, 

i — >C30 

wherever the limit exists. 

Proof. Notice that since each Ui{t) e L'^(M") and ||ui(t)||oo < B, the dominated 
convergence theorem allows for each x E M" 

/ lim u'^{T,x)dT = lim / u'^{T,x)dT 
Jo i^°°Jo 

= Yiva (ui{t,x) — Ui{{),x)) 

i — >oo 

— u{t, x) — u(0, x) for almost every x. 

Hence, by differentiating in i, 

u'(t,x) = lim u[{t^x) for almost every x. 

i — >oc 

□ 



3. The implicit-explicit approximation 

In this section, we consider the case of a 1-dimensional spatial domain, that 
is, a; G M. There is no obstruction to extending any of these results to higher 
dimensions, though it complicates the exposition unnecessarily. 

As is usual, the first task is to define the function spaces to be used. Initial 
conditions will be drawn from a subspace of (M) n L°° (M) , as suggested by Lemma 
[21 and the first four spatial derivatives will be prescribed, for use in Lemma [TUl 

Definition 6. Let 

W = L^(M)nL°°(M)n{/ e C°°(R)|/ has bounded partial derivatives up to fourth or 

For the remainder of this paper, we consider the case where each of the coefficients 
a, e W. Then let X = {/ S W^|5i(||/||i) < 00 and <?oo(|i/|ioo) < 00}. We consider 



the case where the initial condition is drawn from X. 

An approximate solution given by the implicit-explicit iteration will be the piece- 
wise linear interpolation through the iterates computed by ([2]) . A smoother approx- 
imation will prove to be unnecessary, as will be shown in Lemma [TT] 



IMEX METHOD CONVERGENCE FOR A PARABOLIC EQUATION 

Definition 7. Suppose /o and h > are given. Put 

(8) fn+l = {I- hA)-\fn + hG{fn)). 

The function 



(9) 



u{t) =[\-(^-- n{t)Yj fn(t) + - /n(t)+l 



where n{t) = [^J, is called the implicit-explicit iteration of size h beginning 
at /o- 

Calculation 8. We explicitly compute the operator (7 — hA)~^ using Fourier 
transforms. Suppose 

(7 - hA)u{x) = u{x) - hAu{x) = f{x). 
Taking the Fourier transform (with transformed variable w) gives 

u{u>) + huP'u{uj) = f{uj), 



u{co) = 



The Fourier inversion theorem yields 

// 1 p iui{x-y) \ 



Using the method of residues, this can be simplified to give 

(10) u{x) = {{I - hA)-^f) i^) = ^j /(y)e-l^-^l/^dy. 

Calculation 9. Bounds on the and operator norms of (7 — hA)^^ are now 
computed. First, let / e L°°{M.). Then 

\{{I-hA)-'f){x)\ = 



< 

< 
< 



f(y)e-\v-\/^dy 
vh Jo 



so ||(7-/iA)-i||oo < 1. 
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Now, let / e Li(R). So then 

/■° 

||(/~/^A)-V||l = / 

J —i 



2Vh 



< 



1 



oo 

OO /'OO 



dx 



2^/h J- 
1 



\f{y)\e-\y-\/^^dydx 



OO <J —OO 



< I \f{y)\ I e 



< 



which means ||(/ — hA) < 1. 

The thhd condition of Lemma [3] is a control on the slope error of the approxi- 
mation. A bound on this error may be established for the implicit-explicit iteration 
as follows. 

Lemma 10. Suppose fo € X , h > 0. Put f{x, t) ~ foix) + tD{x), where 

{I-hA)-'{fo + hG{fo))-fo 
h 

Then for every < t < h, 

(11) \\f'{t)-F{fm\oo^o{h). 

Proof. Recall every function in X will have bounded partial derivatives up to fourth 
order from Definition [S] 



\f{t)-F{fm[ 



< 



< 



< 



\\D~{A{ff,+tD) + G{fo + tD))\\^ 

COO 
A{fo+tD)+Y,a^ifo+tDy 

D-Afo- tAD - £ a J ^ (') f^ (tD) 



i=0 \j=0 



D-Afo-tAD-J2<^do 



i=0 



Oik) 



(/~/iA)-i -/ 



fo - A/o 



+((/-/iA)-i-/)G(/o)|L +0(/i) 



Now, using the fact that (J - hA)'^ - I = {I - hA)-'^{hA), 
||/'(t)-F(/(0)||oo < ||(/-/iA)-iA/o-A/o 

+ iI-hA)-'ihA)Gifo)\\^ + Oih) 

< \\{I- hA)-\hA){Afo + G(/o))||oo + Oik) 

< /i||(/-/iA)-i(AF(/o))||oo + 0(M 

< h\\{I - hA)-^\U\{AF{fo))\\oc + Oih) ^ 0{h) 
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□ 



Lemma 11. Suppose < ft-^ — > 0. Let Ui be the implicit- explicit iteration of size hi 
beginning at /q £ X on t £ [0, T] . Then provided there exist A,B>Q such that for 
each i and t G [0,T], < A and \\ui(t)\\^ < B, then the sequence 

converges pointwise to a function in t. The limit function is piecewise difjerentiable 
in t. 

Proof. Let Ui be the implicit-explicit iteration of size hi. By Lemma fTOl the slope 
error is bounded: 

\\u[(t) - F{udt))\\oo ^ 0{h,) ^ e,. 
Notice that 0. Then, since X C C^(R"), Lemma [3] applies, giving a pointwise 
limit function u{t). Finally, since the slope error uniformly vanishes. Lemma [5] 
implies that the solution is piecewise differentiable. □ 

4. "yl priori ESTIMATES" FOR THE APPROXIMATE SOLUTIONS 

Now we demonstrate that the implicit-explicit method converges for all initial 
conditions in X. Specifically, for each fo £ X, there exist A,B>Q such that for 
each i and t E [0,r], ||ui(t)||i < A and ||ui(i)||oo < B, given sufficiently small T. 
We begin by recalling that from Calculation [9l the L°°-norm of (/ — hA)^^ is less 
than one. This means that for the implicit-explicit iteration. 



l/n 



-|-l||oo 



< 



< 



< 



< 



< 



fn + hGifn)\\ 



/nil 



E 

i=0 



0-ifn 



fn\\oc. + hJ2hi\\<>c-\\fn\ 
i=0 

oo 

/niloo + ^E ll"«llooll/nl 
fn 1 1 oo hgoo ( 1 1 fn \ \oo ) 



Hence the norm of each step of the implicit-explicit iteration will be controlled by 
the behavior of the recursion 

(12) fn+l = fn + hgooifn), 

for fn, h, a > 0. Since we are only concerned with short-time existence and unique- 
ness, we look specifically at h ^ T/N and < n < iV, for fixed T > and iV G N. 

Remark 12. The recursion defined by (fT^ is an Euler solver for 

(13) ^=.goo(2/), withy(0) = /o. 

This equation is separable, and goo is analytic near /o, so there exists a unique 
solution for the initial value problem p3p for sufficiently short time. Also, whenever 
y{t) > 

d^V , / / 

^ g'ooiyit)) > 0, 

the function y{t) is concave up. As a result, the exact solution to p3p provides an 
upper bound for the recursion (jl2p . More precisely, we have the following result. 
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Lemma 13. Suppose y{0) = /o > m p3|l . Let T > be given so that y is 
continuous on [0,T], and let N CzN. Then for each < n < N , fn < viT), where 
fn satisfies p2p with h = T/N . 

Proof. Since the right side of is strictly positive, the maximum of y is attained 
at T on any interval [0,r] where y is continuous. Furthermore, since y{0) > 0, 
it follows from Remark [T^ that y is concave up on all of [0,T]. Therefore, y is a 
convex function on [0,T]. Hence Euler's method, p^ . will always underestimate 
the true value of y. Another way of stating this is that 

/„ < yinh) < y{T). 

□ 



Using Lemma [131 the growth of iterates to p2)) may be controlled independently 
of the step size. This provides a uniform bound on the sequence of implicit-explicit 
approximations. 

Lemma 14. Suppose < hi = T/i for i G N. Let Ui he the implicit- explicit 
iteration of size hi beginning at fo d X on t G [0, T] . Then there exists a B > 
such that for each i and t € [0,T], we have ||wi(i)||oo < B for sufficiently small 
T > 0. 



Proof. Suppose fi„ is the n-th step of the implicit-explicit iteration of size hi. If 
we let y{0) = ||/o||oo, Lemma [T^ implies that for any i and any < n < i 

WUnWoo <y{T) 

for sufficiently small T. Hence by ([9]) and the triangle inequality, ||ui(t)||oo < B for 
alH and i e [0,T]. □ 

With the bound on the suprema of the approximations, we can obtain a bound 
on the 1-norms. 

Lemma 15. Suppose < hi = T/i for i G N. Let Ui be the implicit- explicit 
iteration of size hi beginning at fo Cz X on t € [0, T]. Then there exists an A > 
such that for each i and t £ [0,r], we have ||7ii(<)||i < A for sufficiently small 
T > 0. 



Proof. First, notice that Lemma fT4l implies that there is a B > such that for each 
i and t e [0, T], we have ||ui(t)||oo < A for sufficiently small T > 0. Again suppose 
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fin is the n-th step of the implicit-explicit iteration of size hi. Then we compute 



II/. 



i,n+l\\l 



< 
< 

< 

< 

< 

< 
< 



fin I 
fin I 

fin I 

fin I 

fin I 

fin I 
fin I 



i + /i»||G(/„0||i 

OQ 

l+h,J2hkf-n 



k=0 



i + h,Y, / \akftn\dx 

i + /i»Eii/-ii-'ii 

^fc II oo II /mil 1 ~t~ 11^0 111 



fc=l 



l + /iiEl|afe||ooS'= ^ +/ii||ao||i 



fe=i 



^ ^ ■^5oo(S) - -^llaolloo ) + /lillaolli 
i(l + /iiC) + /i,||ao||i 



This recurence leads to 



n~l 



11/ 



in 111 



< 



||/olli(l + h,Cr + E (1 + '^'^)' 



< ||/olli(l + /izC)" + /i,|lao|| 



h^C 



< ( ll/o||i + ^l|ao||i) (l + /i.C)"-i||ao||i 

< ( ll/o||i + ^||ao||i) + -^ll«olli 



< 



ll/o||i + ^ll«o||i) (l + ^)^~^ll«olli 



< (ll/o||i + ^l|ao||i)e^^-^||ao||i=A 

Once again, by referring to ([9]) and using the triangle inequality, it follows that 
\\utit)\\i < B for all i and t G [0,r]. □ 

Theorem 16. Suppose < hi — T/i for i £ N. Let Ui he the implicit- explicit 
iteration of size hi beginning at fo £ X on t £ [0,r]. Then, for sufficiently small 
T > 0, the sequence {ui{t)}°l^ converges pointwise to a function in t. The limit 
function is piecewise differentiable in t. 

Proof. This compiles the results of Lemma [11] Lemma [HI and Lemma [151 D 

Remark 17. These proofs can be generalized further to handle all equations of 
the form 



du{t) 
dt 



^ L{u{t)) + G{u), 



where G is as in ([4]). If the operator L satisfies 
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• L : Li(K)nL°°(M)nC°= (K) L°= (R)nC°= (M) is a sectorial linear operator 

• \\{I-hL)-% < 1 and ||(/-/iL)-i||eo < 1, 
then the implicit-explicit iteration 

U+l=^{I-hL)-\f^ + hG{fn)) 

converges for whenever f ^ X. 

5. Conclusions 

The convergence proof for the implicit-explicit method presented here has a 
number of advantages. First of all, like all implicit-explicit methods, each approxi- 
mation to the solution is computed explicitly. As a result, a fully discretized version 
(as is standard in the literature) is easy to program on a computer. Theorem [TB] 
therefore assures the convergence of these fully discrete methods. 

However, since the implicit-explicit method presented here is discretized only in 
time, the convergence proof actually shows the existence of a semigroup of solutions. 
As a result, the convergence proof forms a bridge between the functional-analytic 
viewpoint of differential equations, namely that of semigroups, and the numerical 
methods used to approximate solutions. While the existence and uniqueness of solu- 
tions for ll]) has been known via semigroup methods, the proof provided here gives 
a more elementary explanation of how this occurs. In particular, it approximates 
the semigroup action directly. 

References 

[1] Uri M. Ascher, Steven J. Ruuth, and Brian T. R. Wetton. Implicit-explicit methods for time- 
dependent partial differential equations. SIAM J. Numer. Anal., 32(3):797-823, 1995. 

[2] Michel Crouzeix. Une methode multipas implicite-explicite pour I'approximation des equations 
d'evolution paraboliques. Numer. Math., 35:257-276, 1980. 

[3] Bernold Fiedler and Arnd Scheel. Spatio-temporal dynamics of reaction-diffusion equations. 
In M. Kirkilionis, R. Rannacher, and F. Tomi, editors, Trends in Nonlinear Analysis, pages 
23-152. Springer- Verlag, Heidelberg, 2003. 

[4] Dan Henry. Geometric Theory of Semilinear Parabolic Equations. Springer- Verlag, New York, 
1981. 

[5] John H. Hubbard and Beverly H. West. Differential equations: a dynamical systems approach. 
Springer, 1997. 

[6] Eberhard Zeidler. Nonlinear functional analysis and its applications 11/ A: Linear monotone 
operators. Springer- Verlag, New York, 1990. 



